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The application of the hyperspherical harmonic approach to the case of non-local two-body po¬ 
tentials is described. Given the properties of the hyperspherical harmonic functions, there are no 
difficulties in considering the approach in both coordinate and momentum space. The binding en¬ 
ergies and other ground state properties of A = 3 and 4 nuclei are calculated using the CD Bonn 
2000 and N3LO two-body potentials. The results are shown to be in excellent agreement with 
corresponding ones obtained by other accurate techniques. 


I. INTRODUCTION 

A number of “realistic” nucleon-nucleon (NN) potentials have been determined in recent years hies , that 
reproduce the deuteron binding energy and fit a large set of TVTV scattering data below the pion-production threshold 
with a x 2 /datum ~ 1. The high accuracy achieved in reproducing TV TV observables has spurred renewed interest in 
testing these potential models in systems with A > 3. Several methods have been developed to accurately solve the 
Schrodinger equation for these systems, and meaningful comparisons with precise experimental data, for example for 
the TV — d reaction, are nowadays possible j4 s 5]. Among these methods, of particular importance are the Faddeev- 
Yakubovsky (FY) equations approach 1111 HI fllll . the quantum Monte Carlo methods mm , and the 
techniques based on the exp ansion of the nuclear wave function on an appropriate basis, like the hyperspherical- 
harmonics (HH) [3 (3 OH EH 0, the harmonic oscillator mm, or the gaussian basis 0i- Some of these 
methods are variational while others are not, and they are more or less advantageous depending on the problem at 
hand. All of them provide very accurate results for bound states [23]. On the other hand, the application to scattering 
states presents a number of difficulties that are specific to each of these methods. Of course, scattering wave functions 
are an essential input to the calculation of scattering observables as well as of cross sections for electroweak reactions, 
such as, for example, the low energy radiative and weak capture processes relevant in nuclear astrophysics. The FY 
approach has been developed both in coordinate and momentum space and, in general, gives very accurate results. 
However, this approach becomes problematic when considerin g s cattering processes with charged particles at low 
energies, although recent improvements have been made in Ref. [ll| to reduce some of these problems. In this respect, 
variational methods do not encounter any difficulty, since the electromagnetic interaction between the particles can be 
taken into account in configuration space [3] ■ For example, with the HH technique p — d elastic scattering observables 
have been calculated with the same degree of accuracy as corresponding n — d ones J3 EE B3 ■ 

The failure of theory based the “realistic” TVTV models in predicting successfully A =3 and 4 bound and scattering 
state properties provides evidence for the need of including a three nucleon interaction (TNI) [3H- However, the 
current understanding of the TNI is still in an early stage, and discrepancies between theory and experiments are still 
observed, as in the case of the TV — d [3 EE EE 32] and p — 3 He EE El] A y polarization observable (the so-called 
“A y - puzzle”) It is not clear if these discrepancies can be solved by using more sophisticated models of TNIs, or if the 
problem still resides in the TV TV interaction (in particular, in P waves) I.'IT 13a . 1371] . 

More recently, new potentials have been derived using chiral perturbation theory at increasing order EE El • Also 
many-body forces can be derived on the same footing EE EH EH EE Af- Furthermore, there are other models which 
are constructed up to a certain cutoff momentum A ( low-q potentials | l5l l4fi| ). All these new potentials are defined 
primarily in momentum space and are non-local. Therefore, accurate techniques which can solve the corresponding 
Schrodinger equation for A =3 and 4 are important. At present, only the FY I El and the No Core Shell Model |20l| 
methods have been employed for this task. In this paper, we will show that the HH method too can be successfully 
applied to treat this kind of potentials. 

In the last few years, considerable effort has been devoted by the authors of the present paper to the development 
and application of the HH technique to study bound and scattering states of three or four nucleons with realistic 
local TVTV potentials 114 IlH fl6l 1471 148L l4fl. A version of the method, that has been rather extensively exploited in 
these calculations, includes an appropriate correlation factor in the HH basis functions, so as to take into account the 
strong short-range repulsion of the TVTV interaction and, therefore, improve the rate of convergence of the expansion. 
This approach is known as the correlated-hyperspherical-harmonics (CHH) method EH EH El- Up until now, the 
calculations have been done only for local nuclear interactions, in particular the Argonne uis two-nucleon (AV18) 0] 
together with the Urbana IX (UIX) three-nucleon El interaction. 
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The HH method, without correlation factors, has been also employed for the calculation of trinucleon 0 and a 
particle 0 bound states and four-nucleon scattering processes |5Cj]. Such a method appears especially convenient in 
the case of non-local potentials. The object of the present paper is the calculation of the trinucleon and a particle 
bound states with the uncorrelated HH expansion using two recent, realistic non-local two-nucleon interactions, 
namely the CD Bonn 2000 3 and the N3LO Ohl potentials. This kind of expansion, in fact, can be performed 
equally well in coordinate or momentum space j.5l| . It is thus possible to treat in the corresponding space the part 
of interaction given in coordinate or momentum space. With two-body local potentials the matrix elements of the 
interactions can be obtained via two dimensional integrals. The non-locality of the interaction merely requires a 
three-dimensional integration, which can be performed using standard numerical methods. Obviously, there are no 
difficulties in including the local Coulomb interaction. The same is valid for the most commonly used TNIs, which 
are local interactions. However, the TNI has not been considered in the calculations reported here. As a final remark, 
it should be noted that the present work is only the “starting point” for the implementation of the HH method in the 
study of A =3 and 4 scattering states with non-local interactions. 

The paper is organized as follows. The details of the formalism are presented in the next section. In Sec. IIIII the 
results for the binding energies and ground state properties of systems with A =3 and 4 are presented and compared 
with those obtained with other different techniques. A few final remarks and conclusions are given in Sec. H3 


II. FORMALISM 

In general, an A-body bound state can be written as 

|*(l,2,--.,A)>=£c, 1 |* #t ), (2.1) 

where I'&p) are a suitable complete set of states, and /i is an index denoting the set of quantum numbers necessary to 
completely determine the basis elements. The coefficients of the expansion can be calculated using the Rayleigh-Ritz 
variational principle, which states that 


(6 C *(1,2,---,A)\H-E\*(1,2,---,A))=0 , (2.2) 

where 5 C \&(1, 2, • • •, A) indicates the variation of (1, 2, • • -, A) for arbitrary infinitesimal changes of the linear coeffi¬ 
cients c^. The problem of determining c p and the energy E is then reduced to a generalized eigenvalue problem, 

^<^|J?-£7|^)V=0. (2.3) 

A 

The main difficulty of the method is to compute the matrix elements of the Hamiltonian H with respect to the basis 
states |4' /J ). Usually H is given as a sum of terms (kinetic energy, two-body potential, etc.). The calculation of the 
matrix elements of some parts of H can be more conveniently performed in coordinate space, while for other parts 
it could be easier to work in momentum space. Therefore, it is important that the basis states |'I , M ) have simple 
expressions in both spaces. The HH functions do have such a property, as will be shown below. 

Let us first consider the expression of the HH functions in coordinate space. The internal dynamics of a system of 
A identical nucleons of mass m is conveniently described in terms of a set of N = A — 1 Jacobi vectors X\ pi ..., Xn p , 
constructed from a given particle permutation denoted with p, which specifies the particle order i,j, k ,.... In par¬ 
ticular, xn p = (rj — n)/y/2 (r.j is the position of the i-th particle, etc.), and p = 1 is chosen to correspond to the 
particle order 1,2,3,.... The hyperradial coordinates are defined to be 


P 

^ P) 


N 


x 2 
^ip 




X2 p, , Xj\j pi (f)2pi i &Np\ 


(2.4) 

(2.5) 


where p is the so-called hyperradius, a set of angular-hyperangular coordinates and the suffix (p) recalls the 
use of the coordinate space. Note that p does not depend on the particle permutation used to construct the Jacobi 
vectors. The angles (j> 2 P , • • •, (j)Np are the hyperangles, defined as [53] 


tctn (pip — 


°xp 


\ 


Z. x jv 
J=1 


i = 2,...,N . 


( 2 . 6 ) 
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In terms of these variables, chosen any particle permutation p , the kinetic energy operator becomes 

T = V V 2 ™-l d _ A 2 (fIp P) ) 

2 m .~^ N x<1 ’ 2 m y dp 2 p dp p 2 

The HH functions Y[g](^p^) are the eigenfunctions of the operator A 2 (flp P ' 1 ), and their explicit expression for a generic 
A-nucleon system can be found, for example, in Refs. [5lli£2l Here we consider in some details only the A=3 case. 
The Jacobi coordinates for A= 3 nucleons are defined as 

x 2 P = -^(rj-n) , 

x ip = Y f ( r fc ~ r i + r i)) > ( 2 - 8 ) 

where p = 1 corresponds to the order 1,2,3. The HH function (ilp > ' > ). with a definite value of the total orbital 
angular momentum L, L z , can be written as 0 : 



Y^(n^) = 


( X 2p ) 0 ( X 1 p) 


N[g] (cos (f> P Y 2 (sin (ftp)* 1 P, 


n^+Iib+| 


(cos 2 cf)p 


(2.9) 


J LL Z 


Note that in the A = 3 system only the hyperangle fa P is present (indicated with cp p ). In Eq. <12.911 . Yt^xip) and 
Y( 2 (x 2 p) are spherical harmonics on the two internal Jacobi coordinates Xi p and X 2 V , N[g] is a normalization factor 

and P n2 2 ’ 2 (cos 2 (j) p ) is a Jacobi polynomial, ri 2 being the degree of the polynomial. The grand angular quantum 

number G is defined as G = 2 n 2 + t\ + H 2 . The notation [G] stands for [l\, f? 2 ; n 2 ]. 

Moreover, we consider the anti-symmetrized functions 3 ^{g}(^ ( ' p ^) given by the product of an HH and a spin-isospin 
function. The antisymmetry is obtained by writing J^g}^^) as a sum of terms constructed by starting from all 
the possible even permutations p of the particles, assuming antisymmetry in the interacting pair. For A = 3 they are 
explicitly given by 


y{G}(n {p) ) = 

v 



( 2 . 10 ) 


where the spins (isospins) of particle i and j are coupled to S 2 (T 2 ), which is itself coupled to the spin (isospin) of 
the third particle to give the state with total spin S (isospin T, T z ). The total orbital angular momentum L and the 
total spin S are coupled to the total angular momentum J, J z . Each set of quantum numbers {G, G, L, S 2 , T 2 , S, T} 
is called “channel”, and here the notation {G} stands for 


{G} = {£ 1 J 2 ,L,S 2 ,T 2 ,S,T-n 2 } . (2.11) 

The requirement of the antisymmetry of A’{g}(^^ p ' ) ) constraints the possible choices of the set {G} to those for which 

(_)*2+S 2 +t 2 = _ (2.12) 

The corresponding expressions for the A = 4 system can be found in Ref. 112 - 

In this paper, we will consider modern two-body potential models which act on specific spin and angular momentum 
states of the two-body system. Due to the presence of the sum over the permutations in the expression for 3 Ag}> 
a given particle pair is not in a definite angular and spin state. However, the HH functions with the grand angular 
quantum number G constructed in terms of a given set of Jacobi vectors X\ p , ■ ■ ■ ,xn p , defined starting from the order 
of particles, can be always expressed in terms of the HH functions constructed, for instance, in terms of 
®i( p =i), ■ • •, x N (. p= i ). For example, when A = 3, the following relation holds 


\rLL z 
1 [£i ,(- 2 \n 2 ] 



E 


/~ L L z 

^1,^2,^2 ; [GJ 2 ’ n 2 ] 


(H 


(p) \ 

( P =i)J > 


(2.13) 


where the HH functions are defined in Eq. iTT9l) and the sum is restricted to the values £[, i' 2 , and n 2 such that 
£\ +1' 2 + 2n 2 = G. The coefficients af^' e t , ^ n , relating the two sets of HH functions are known as the Raynal- 
Revai coefficients [53j. There exist several procedures to compute these coefficients for A = 3 and 4 systems |5^[ 13 
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mm min. For A = 3, they can be computed rather easily using the orthonormality property of the HH functions, 
namely 


a 


(p), L 


/ 


dfl 


(p) 

(p=i) 


( 


~yLL z 
1 \tiM ;n 2 ] 



\rLL z 



(2.14) 


For A = 4, we have used the procedure devised in Ref. 0], where the corresponding coefficients are obtained by using 
a set of recurrence relations. Also the spin-isospin states can be recoupled to obtain states where the spin/isospin are 
coupled in a given order of the particles. The result is that the antisymmetric functions Wg} can be expressed as a 
superposition of functions constructed in terms of a given order of particles i, j, fc, . . ., each one having the pair i.j in 
a definite spin and angular momentum state. When the two-body potential acts on the pair of particles i,j, the effect 
of the projection is easily taken into account. 

We now consider the expansion states I’Fp) of Eq. ED- In coordinate space, they have been chosen to be given by 


= fi(p)y {G} (rt p) ), ( 2 . 15 ) 

where p stands for the set of quantum numbers ({G},/), and fi(p) for l = 1,... is a complete set of hyperradial 
functions. Note that the index p has been suppressed, since is at this stage independent on the permutation. 
The corresponding states in momentum space can be obtained as follows. Let hki p , • • •, hk^ p be the conjugate Jacobi 
momenta of the Jacobi vectors. Explicitly, in the A=3 case, these momenta are 


hk‘ 2p 

hki p 



Pi) , 



; iPi+Pj )) > 


Pi being the momentum of the i -th particle. The following relation then may be shown to hold |5lll52l |: 


(2.16) 


( p? • • • i *^Np | • • • i kNp ) — 

= (2„)3A72 1 A C." ' x jp) = 

where C = G + (3 N — 3)/2, Jc+i/ 2 {Qp) is a Bessel function, il p > ' 1 is given in Eq. (12.511 . and 


Q = 


N 



i =1 


(2.17) 


(2.18) 


is the hypermomentum, which plays the role in momentum space corresponding to p in coordinate space. Also Q 
does not depend on the particular permutation used to construct the vectors hk\ p , • • •, hk^p- The momentum-space 
angular-hyperangular variables are defined as 

— [^ip, k^p , , k^p, P2pi > PNp] ; (2.19) 


and 


tcLIl (flip — 


kr 


ip 


\ 3=1 


Then, the momentum space version of the wave function given in Eq. (12.1511 i 


is 


( 2 . 20 ) 


(fci, • • • ,k N | 


/ 


dxi p • 


dxjsip 


p kjp’ x jp 

p^m-2 /<(»)r IG| (nM). 


( 2 . 21 ) 
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Note that also ( k \,..., k jy | ) cannot depend on the permutation index p. In fact, using the expansion of Eq. (I2T7I) 

and the orthogonality property of the HH functions, it follows that 

fcAr | ) = g G AQ)y { G}(rt Q) ) , (2.22) 

where is the same as (V{G}(^^) with x ip —> k ip , and 

(■oo 3AT_i 

gcAQ) = H) G J o d P (q^ 3 jv/2-i J c+i/2{Qp) flip ) • (2.23) 

In the present work, we have adopted two choices for the functions flip), for which the Fourier transform given in 
Eq. 12.2311 can be obtained (almost) analytically. The two choices are illustrated below. 

1. Exponential functions of the type 


flip) = e~ aiP , (2.24) 

where a\ = (3 and cq = (3 + j 1 ^ 1 for l = 2, • • • M, (3 and 7 being two non-linear parameters, which need to be 
optimized depending on the used nuclear potential model. As an example, when A = 3, /? = 0.78 fm -1 and 
7 = 1.2 fm ” 1 for the CD Bonn 2000 potential. Other choices of cq have been made, for instance cq = (3 + J (Z — 1) 
or cq = /?7 _1 , but no significant differences have been found in the final results, once the non-linear parameters 
/3 and 7 are optimized and convergence on M has been reached. The advantage of using the exponential 
basis is that the corresponding functions g G ,i(Q ) have an easy expression, and are given by: 


(n) - ( ,-)G T i D + G ) 

9G,iiQ) = i~ l ) 


,d 


, ,D 


P; 


2 \-£-± D/2 

(1 — U z ) * 2 


1 —G —-D/2 


A) , 


1 + Q 2 /a 2 


(2.25) 


where D = 3 N and P™ is an associated Legendre function. For A = 3 this expression can be written in terms 
of analytical functions, namely 


9G,i(Q ) 


(-i) G \/5! d 3 + Q 2 - cq) G+2 " 

QG+ 1 JW dtf [ y/a* + Q' 2 . 


(2.26) 


2. Another useful form for flip), adopted for example also in Ref. is 

(jD_ 1 ) (7 P) e_ ^ P 


flip) =l D/2 \ 


l\ 


L 


H + D-l)\ 1 


(2.27) 


where d!) D 1 \ i yp) are Laguerre polynomials. Here, there is only one non-linear parameter, 7 , to be variationally 
optimized. In particular, 7 can be chosen in the interval 3.5-4.5 fm -1 for the CD Bonn 2000 potential and 6-8 
fm -1 for the N3LO potential, for both A = 3 and 4. The corresponding functions g G ,iiQ) are less trivial to 
calculate, and are given by: 


**(«> = £>* 2 ‘ +D r < G + fc + D > (1 CT/AV) 


(2.28) 


where u = 


\ A +(2 q / 7 ) = 


and b l k are given by 


bl = 


(-l) fc fl + D-1 


k\ 


l - k 


(2.29) 


so that L\ \x) = foj xk H3- Studies on the convergence on M for the A=3 and 4 calculations are 

presented in the Sec. Hill 


In summary, the variational state is given by 


M 

I ^( 1 , . . . , A) ) = y,y>G>, ; | : 

{G} 1=1 


(2.30) 
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where, in momentum space the expansion states are 

(fei,... ,fejv | ^{G},*) = 9G,i{Q)y{G}(^ Q ^) t (2-31) 

and in coordinate space are 

(x ll ...,x N \y {G} , l ) = M P )y {G} (n ^). ( 2 . 32 ) 

These two expressions can be used to evaluate the matrix elements in Eq. iTOl) . In particular, the normalization 
(TV) and kinetic energy (T) operator matrix elements can be computed both in coordinate and in momentum space. 
Explicitly: 


N {G'},l'-,{G},l 


J{G'},Z';{G},I 


dpp 3N_1 fi'(p) fi(p) J dn^y* {G , } (n^)y {G} (n^) 

dQQ 3N - 1 g G ',i'(Q)9G,i(Q) f dn^y* {G , } (n^)y {G} (n^) 


-Lf 


r d 2 3N — Id G(G + 3N-2) 

dp 2 ^ p dp p 2 


fi(p) 


x j dn^y* {GI} (n^)y {G} (n^) 
J dQQ 3N+1 g G ,^(Q)g Gtl {Q) 

x f dn^ y* {G , } (n^)y {G} (n^) 


(2.33) 


(2.34) 


where it has been used the fact that the HH functions are eigenfunctions of the operator A 2 (flp P ^) defined in Eq. ra 
corresponding to the eigenvalues G(G+3IV—2). In Eqs. (12.3311 and <E3U), as well as in the rest of the present work, the 
permutation index p is omitted, and the integration variables corresponding top = 1 are used (i.e. ^j). 

The calculation of the two-body potential energy matrix elements is more conveniently performed in either coordi¬ 
nate or momentum space, depending on the particular potential model of interest. First of all, due to the antisymmetry 
of the wave function, the following relation holds 


V{G'},l'-,{G},l = ('^{G'},Z'l^l’I'{G},z) 


A(A-l) 

2 


('J , {G'},Z'Kl: 2 )l'I , {G},z) , 


where u(l, 2) acts on the particle pair 1,2 and can be written as 


u(l,2) = V(x' n -,x n ) , 


in coordinate space, and 

^(1)2) = V(k' N \kjy) , 


(2.35) 


(2.36) 


(2.37) 


in momentum space (clearly V and V are related by a Fourier transform). Then, 


V{ G '},i'-{ G },i = AiyA 2 ^ J dx 1 ---dx N J dx' N fv{p l )y* {GI} {^ p ' ) ) 

x V(x' N ; x N )fi{p) y {G} (fl (p) ) 


(2.38) 


or 

V {G >},1';{G},1 = A(A ~ j dk,---dk N J dk' N g G , tl '(Q')y* {G , } (nW) 

xV(k' N -, k N )g G ,i(Q) ^{G}(^ (Q) ) , (2.39) 

where p’ , ) are the hyperradial coordinates associated to the Jacobi vectors Xi, ■ ■ ■, x^ N _ 1 ),x' N , etc. When the 

potential energy operator is a local operator in coordinate space, namely 

V(x' n ]x n ) -> Vioc( x n )6(x' n - x N ) , 


(2.40) 
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it is more convenient to calculate the corresponding matrix elements using Eq. which simplifies to 


J dxi -■ ■ dx N fv{p)y* {G , } (n < ' p) )Vi oc (x N )fi{p)y {G }{^ p) ) . ( 2 . 41 ) 

Examples of two-nucleon potential models of this form are the Argonne vis [lj or the Nij megen II pj potentials. Note 
that previous applications of the HH method were limited to these cases in Refs. jjjj,[l7|. On the other hand, for the 
CD Bonn 2000 or N3LO potentials, which are non-local operators in momentum space, the use of the momentum- 
space expression of V{g'},z';{G},z is more convenient. In general, the potential energy operator is given as a sum 
of different local and/or non-local terms. The computation of the matrix elements of each part can be performed 
using either the coordinate- or the momentum-space expression, depending on the convenience. For instance, if the 
selected potential energy model includes the CD Bonn 2000 potential and the (point) Coulomb interaction, the matrix 
elements of the CD Bonn 2000 are computed using Eq. (J239, the ones of the Coulomb potential are computed u sing 
Eq. 12.4111 . The same procedure can be applied to the TNIs mostly used in the literature, i.e. the Urbana-type [lj] 
and Tucson-Melbourne-type 0, 0] potentials. However, since the aim of this work is to study the applicability of 
the HH method when non-local potentials are used, the TNI has not been included. 

Other potentials models which are frequently used in the literature are those developed by Doleschall et al. 1110 ]. 
These models consist of non-local operators given in coordinate space. In that case, it is more convenient to perform 
the calculation using Eq. 12311). _ _ 

The calculation of the integrals involved in Eqs. 123E1> or 1239 is a non-trivial numerical task. As an example, 
let us consider Eq. 1239- Remembering that V(k' N ; kpj) acts on the particle pair 1,2, it is convenient, using the 
Raynal-Revai coefficients and the recoupling of spin-isospin states, to express the states (4 '{g , },z , I and l^'{G},z) m 
momentum space as a superposition of HH functions and spin-isospin states constructed using the order of particles 
1,2,3,..., A. Most of the integrations can be now performed analytically and the matrix elements, for a general 
A-body system, reduces to a sum of three-dimensional integrals of the type 


1 G' X N , n ' N \G,f.N,n N 


where 


/• OO /* OO /*00 

/ dqq D ~ 4 dk N (k N ) 2 / dk' N (k' N ) 2 g G ',v(Q') 
/ 0 Jo Jo 


x(cos <p’ N Y' N (sin <p' N ) v 'P”i+ D/2 5/2 ' eN+1/2 {cos2ip , N )v 1 e , NtS ,. eNtS (k' N ,k N ) 
x 9G,i(Q) (cos <PnY n (sin p N y P^ D,2 ~ b,2/N+1,2 (cos 2 ip N )S v y , 


(2.42) 


Q 2 = k% + q 2 , ( Q') 2 = ( k' N ) 2 + q 2 , cos ip N = k N /Q , cos ip' N = k' N /Q' , 

and 

v = G — 2riN — (-n i v' = G' — 2n' N — i' N . 


(2.43) 

(2.44) 


Moreover, in Eq. *E 22 J) v r N ,s'-,e N ,s(k'N’ k]y) is the two-body potential acting between states of the pair of particles 
1,2 of total angular momentum j, and orbital angular momentum and spin quantum numbers P N ,S' (on the left) 
and £n, S (on the right). The integrals IcYi' ^-nZ-G e N n N can computed beforehand and stored in computer disks. 
The last step consists in combing dcYf/^nfc t N n N w ith the Raynal-Revai coefficients to obtain the matrix elements 
of Eq. (12.391) . Finally, the inte grat ions involved in Eq. i|2321) can be accurately performed with standard numerical 
techniques (Gauss integration) |60j . 


III. RESULTS 

In this section, the binding energy and ground-state properties obtained for the nuclear systems with A=3 and 4 
are presented for the CD Bonn 2000 3] and the N3LO [3E| momentum space potentials. These interaction models 
are decomposed on partial waves, and the partial wave decomposition is truncated at a certain value of the two-body 
total angular momentum j m ax ■ In the present work j ma x = 6 has been chosen, which allows for an accuracy of better 
than 1 keV in the triton binding energy. 

The section is divided into three subsections: in Sec. IIII Al and IIII Bl the convergence of the expansion with respect 
to the quantum numbers {G} and the number M of hyperradial functions is discussed for A =3 and 4, respectively. 
The converged results for the ground-state properties of triton, 3 He and 4 He are presented in Sec. mm and compared 
with the results obtained with different approaches. 






TABLE I: Triton binding energies in MeV, calculated with the CD Bonn 2000 two-nucleon interaction, using only the first 
3 expansion channels and the exponential functions of Eqs. i pm and <2.2(11 as expansion basis for the hyperradial and 
hypermomentum functions. M is the maximum number of basis elements, G is the grand angular momentum. 


G 

00 

II 

£ 

M=10 

M= 12 

r—1 

II 

§ 

M=16 

40 

7.6764 

7.6835 

7.6842 

7.6845 

7.6846 

50 

7.6845 

7.6921 

7.6934 

7.6934 

7.6935 

60 

7.6870 

7.6948 

7.6958 

7.6962 

7.6963 

70 

7.6880 

7.6958 

7.6969 

7.6974 

7.6976 

80 

7.6885 

7.6964 

7.6974 

7.6979 

7.6981 


TABLE II: Same as Table [I] but with the Laguerre polynomial expansion of Eqs. STT7\ and i THSI for the hyperradial and 
hypermomentum functions. 


G 

40 

50 

60 

70 

80 


M=8 M=10 M= 12 M=14 M=16 M=18 A/=20 M=22 M =24 


7.6582 7.6773 7.6827 
7.6665 7.6862 7.6916 
7.6691 7.6892 7.6946 
7.6702 7.6904 7.6959 
7.6706 7.6910 7.6965 


7.6835 7.6842 7.6846 
7.6924 7.6931 7.6936 
7.6955 7.6962 7.6966 
7.6967 7.6974 7.6978 
7.6973 7.6980 7.6984 


7.6847 7.6848 7.6848 
7.6937 7.6937 7.6937 
7.6967 7.6967 7.6967 
7.6979 7.6980 7.6980 
7.6985 7.6986 7.6986 


A. Convergence of the A= 3 Results 

In this subsection, we show the level of accuracy reached by the method presented in this work. As an example, in 
Table III and IT11 the triton binding energy is calculated with the CD Bonn 2000 potential using only the first 3 channels 
of Table I of Ref. 0 , and increasing the value of the grand angular momentum G and the value of the number M of 
basis elements in the expansion of the hyperradial and hypermomentum functions fi(p) and gG,i(Q ) (see Eqs. (12 .1 511 
and 12.2211 L In Table in the exponential basis of Eqs. m> and m> is used, and the non-linear parameters (3 and 
7 are 0.78 fm -1 and 1.2 fm -1 respectively. In Table gh up) and gc,i(Q) are expanded on Laguerre polynomials, 
Eqs. (l2~271) and and the non-linear parameter 7 is 4.0 fm 1 . 

By inspection of the tables, we can conclude that i) the convergence of the exponential basis at the 1 keV accuracy 
level is slightly faster than the one of the Laguerre polynomials basis. In fact, the maximum value of M used in the two 
expansions, M max , is 16 and 24 respectively. However, it should be noted that the exponential basis gives numerical 
problems for higher values of M max , and turns out to be not suitable for the N3LO potential. Nevertheless, the 
binding energy B at G=80 calculated with M max = 16 exponentials is only 0.5 keV smaller than the result obtained 
with M max = 24 Laguerre polynomials, ii) In the case of the exponential basis, the expansion on M is truncated when 
B{M max ) — B(M max — 2) < 0.5 keV. In the case of the Laguerre polynomial expansion, BiM max ) — B{M max — 2) 
is smaller than a tenth of keV. iii) The convergence on G is reached at G max = 80. In fact, in Table QJ we can see 
that going from G=60 to G=70 and from G=70 to G=80, the gain in B is 1.3 keV and 0.5 keV, respectively. In 
Table El the corresponding values are 1.3 keV and 0.6 keV, respectively. Due to these considerations, we are confident 
in concluding that the results of 7.699 MeV for the triton binding energy, obtained with only the first 3 channels in 
our channel expansion, is accurate at least at the 1 keV level. Note that the value for G max used in Ref. [0 for the 
AV18 potential and the first 3 expansion channels is 180. This is related to the fact that both the CD Bonn 2000 and 
N3LO potential models have a rather soft short-range repulsion compared with the AV18 potential model. 

A similar procedure to reach the convergence on G and M has been used also when other channels are included in 
the calculation. The converged A= 3 results presented in Sec. irrm have been obtained including up to 23 channels, 
18 having total isospin T= 1/2 and 5 having T= 3/2 (see Sec. Inland Ref. HH). In particular, for the CD Bonn 2000 
(N3LO) potential, the maximum value of the grand angular momentum G max is 80 (40) for the first 3 channels, 40 
(20) for the next 5 ones, and 30 (20) for the last 10 channels with T=l/2. The convergence for the T= 3/2 channels 
is reached using G rnax = 20 . 


B. Convergence of the A=4 Results 

First of all, we discuss the convergence of the expansion with respect to the quantum numbers {G}. We follow 
the procedure adopted in Ref. 0 The states ^{g},;) are separated in six classes , depending on the total isospin 










9 


TABLE III: Convergence of a-particle binding energy B (MeV), mean value of the kinetic energy (T) (MeV), mean square 
radius \J ( r ) 2 (fm) and P— and D— wave percentages for different values of M, the number of functions gG,i{Q) included 
in the expansion. In this example, gG,i(Q ) are obtained expanding the hyperradial functions on Laguerre polynomials with 
7 = 6.0 fm -1 . The potential model considered is the N3LO. The basis includes HH states with {Gi, G 2 , G 3 , G 4 , G 5 , Gg} = 
{24, 20,16,16,16, 20} (see the text for more details). 


M 

B 

(T) 

\/W) 

Pp 

Pd 

8 

24.649 

72.326 

1.453 

0.172 

9.272 

10 

25.195 

69.940 

1.495 

0.173 

9.315 

12 

25.339 

69.309 

1.512 

0.172 

9.287 

14 

25.362 

69.244 

1.514 

0.172 

9.290 

16 

25.373 

69.230 

1.516 

0.172 

9.289 

18 

25.376 

69.236 

1.516 

0.172 

9.289 


quantum number T and on £ 1 , £ 2 , £ 3 , the values of the orbital angular momentum quantum numbers associated to the 
three Jacobi vectors (or momenta). The most important classes are those ones with T = 0 and with the lowest values 
of f'i+^2+^3- The convergence of the class i is studied by including in the expansion all the states | ) with grand 

angular quantum number G < Gi, and then increasing the value of G,;. The rate of convergence depends primarily 
on the repulsion at short interparticle distance of the adopted potential models. More details on this argument and 
the choice of the basis can be found in Ref. Ca¬ 
in the case of the two potential models considered in this work, it turns out that both the models have a rather soft 
repulsion and the convergence is reached for values of Gi, i = 1,..., 6 rather smaller than those found in Ref. na, 
where the the AV18 d and Nijmegen II |2) potential models were considered. To give an example, for those models the 
convergence was reached for {Gi, G2, G3, G4, G5, Gq} = {72,40, 34,28, 24, 20}, while for the CD Bonn 2000 potential 
a satisfactory convergence is obtained for {Gi, G2, G3, G4, G5, Gq} = {52,34,28,24,20,20}. The N3LO potential is 
derived from an effective field theory by means of an expansion valid at low momenta, and it is strongly suppressed 
at high momenta. This corresponds in coordinate space to a rather soft repulsion at short interparticle distances. 
Correspondingly, a satisfactory convergence is reached at {Gi, G 2 , G 3 , G4, G 5 , G 6 } = {24,20,16,16,16,20}. The 
“missing binding energy” due to the truncation of the basis has been estimated, using the same procedure adopted 
in Ref. [lj|, to be around 10 keV for the CD Bonn 2000 case, and around 5 keV for the N3LO potential. 

We now consider the convergence with respect to the number M of functions gG,i{Q) used in Eq. (12.301) . In 
Table EH the binding energy and other ground state properties calculated for the N3LO potential for different 
values of M are reported. In this example, the hyperradial functions have been expanded on Laguerre polynomials 
with 7 = 6.0 fm -1 . All the HH functions belonging to the six classes discussed above, namely for the choice 
{Gi, G2, G3, G4, G5, Ge} = {24,20,16,16,16,20}, have been included. From the table it is seen that all quantities 
have converged quite well for M « 16 — 18. 


C. Results for A=3 and 4 

The converged values for the triton, 3 He and 4 He binding energies and other ground-state properties are listed 
in Table EH El and ED respectively. All these results have been obtained expanding the hyperradial functions on 
Laguerre polynomials, with 7=4 fm -1 and 7 fm -1 for the CD Bonn 2000 and N3LO potentials, respectively. The 
numerical uncertainty in the triton and 3 He binding energy has been estimated to be at the most of the order of 1 
keV, and in the 4 He binding energy of the order of 10 keV. These results are compared with those obtained with 
other approaches. In particular, we have considered the Faddeev (for A= 3) and Faddeev-Yakubovsky (FY) (for A= 4) 
approach 0, HD ED ED > and the No Core Shell Model (NCSM) approach ED ED- For sake of comparison, the 
results for the AV18 potential are also listed pi HU ITtL ED 1691. Note that the HH results of Ref. HD do not include 
the T = 3/2 states. In the 3 He and 4 He case, all the results for the CD Bonn 2000 and N3LO potentials have been 
obtained including the point Coulomb interaction, except for the 4 He FY results with the N3LO potential. The latter 
have been obtained including a more complicated electro-magnetic interaction between nucleons j6p . whose effects 
are however quite similar to those of the point Coulomb interaction. 

There is a good agreement, at the level of 0.1 %, among the results obtained with different techniques for all 
quantities considered in the tables. Furthermore, the mean value of the kinetic energy for the N3LO case is smaller 
than that found with the AV18 potential. This is due to the fact that the repulsion at short interparticle distances is 
softer for the N3LO potential than for the AV18 potential. Also the percentages of the P- and D-waves is significantly 
smaller for the N3LO potential. This explains the faster convergence of the HH expansion in this case. The CD Bonn 
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TABLE IV: The triton binding energies B (MeV), the mean square radii { r 2 ) (fm), the expectation values of the kinetic energy 
operator (T) (MeV), and the mixed-symmetry S', P, D, and isospin T= 3/2 probabilities (all in %), calculated with the CD 
Bonn 2000 and N3LO potentials, are compared with the results obtained within the Faddeev equations approach [Til. Ifijj . Ififl El 
(FE) and within the No Core Shell Model (NCSM) approach l2fl IfiSH . The results obtained in Refs. fiTf and |65] within the 
HH, correlated-hyperspherical-harmonics (CHH) and FE approaches for the AV18 potential have been also reported for sake 
of comparison. These last HH results do not include the T= 3/2 states. 


Interaction 

Method 

B 

(T) 

VW) 

P.S 1 

Pp 

Pd 

Pt= 3/2 

CD Bonn 2000 

HH (this work) 

7.998 

37.630 

1.721 

1.31 

0.047 

7.02 

0.0049 


FE [65, 66) 

7.997 

37.620 

- 

1.31 

0.047 

7.02 

0.0048 


FE [10, .671 

7.998 

37.627 

- 

1.31 

0.047 

7.02 

0.0048 


NCSM [20, 68] 

7.99(1) 

- 

- 

- 

- 

- 

- 

N3LO 

HH (this work) 

7.854 

34.555 

1.758 

1.36 

0.037 

6.31 

0.0009 


FE [65, 66] 

7.854 

34.546 

- 

1.37 

0.037 

6.32 

0.0009 


FE [10, 671 

7.854 

34.547 

- 

1.37 

0.037 

6.32 

0.0009 


NCSM [20, 68] 

7.85(1) 

- 

- 

- 

- 

- 

- 

AV18 

CHH [65] 

7.624 

46.727 

- 

1.293 

0.066 

8.510 

0.0025 


HH ]16] 

7.618 

46.707 

1.770 

1.294 

0.066 

8.511 

- 


FE [65] 

7.621 

46.73 

- 

1.291 

0.066 

8.510 

0.0025 

TABLE V: Same as Tabic EH but for 3 Hc. 

Interaction 

Method 

B 

(T> 

v'F) 

Ps' 

Pp 

Pd 

II 

03 

to 

CD Bonn 2000 

HH (this work) 

7.262 

36.777 

1.759 

1.54 

0.046 

7.00 

0.0109 


FE [65, 66] 

7.261 

36.756 

- 

1.54 

0.046 

7.00 

0.0110 


FE [10, 67] 

7.263 

36.761 

- 

1.54 

0.046 

7.00 

0.0110 

N3LO 

HH (this work) 

7.128 

33.789 

1.797 

1.61 

0.037 

6.31 

0.0062 


FE [65, 66] 

7.128 

33.775 

- 

1.61 

0.037 

6.31 

0.0063 


FE [IQ, £2] 

7.128 

33.775 

- 

1.61 

0.037 

6.32 

0.0063 

AV18 

CHH [££| 

6.925 

45.685 

- 

1.530 

0.065 

8.467 

0.0080 


FE [65] 

6.923 

45.68 

- 

1.524 

0.065 

8.466 

0.0081 


2000 is intermediate between the two cases discussed above. 

We now consider the results for the T=3/2 and T > 0 components in the A=2> and 4 systems, respectively. The 
percentage of the T = 1 component in the 4 Heground state wave function is almost independent from the adopted 
potential model. In fact, as discussed in Ref. [[3, 50% of it is due to the effect of the Coulomb interaction between 
the protons, the remaining 50% is due to the charge symmetry breaking (CSB) terms in the nuclear interaction. 
Consequently, the difference in Pt=i for the various two-body potentials is small. On the contrary, the T=3/2 
component in triton and the T = 2 component in 4 He are largely dominated by CSB of nuclear origin (different pion 
masses, etc.). The values reported in Tabic Hvl and I VII show that, depending on the interaction, rather different values 
for the triton Pt— 3/2 and the 4 He Pt =2 are obtained (note that the standard models of TNI have little effect on the 
isospin admixtures [TtJ]). The CSB magnitude in the two-body interaction is fixed by fitting the NN scattering data, 
and therefore comes from differences observed in the pp and np systems. In particular, the CD Bonn 2000 and N3LO 
potentials fit the same NN data set [Tcj- The origin of the rather large differences found for the triton Pt— 3/2 and the 
4 He Pt -2 (a factor 5 between CD Bonn 2000 and N3LO) must be related to quite different off-shell behavior of the 
CSB terms of both interactions. Note that the T= 3/2 component in 3 He is also affected by the Coulomb interaction, 
which reduces by more than a factor of 2 the difference between the CD Bonn 2000 and N3LO Pt- 3/2 results. 

Finally, we observe that knowledge of the T = 1 and 2 percentages could be important for parity-violating electron 
scattering experiments on 4 He, aimed at studying admixture of strange quark ss pairs in nucleons and nuclei El El 
E3- It could play an important role also in the study of the reaction d + d —> a+ tt°. This reaction is possible only if 
isospin symmetry is violated, namely it probes directly the CSB terms in the nuclear Hamiltonian 1741751 . 
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TABLE VI: The a-particle binding energies B (MeV), the mean square radii yj(r 2 ) (fm), the expectation values of the kinetic 
energy operator (T) (MeV), and the P, D, T = 1 and T = 2 probabilities (%) for the two non-local potentials considered in 
this paper. The results obtained in Ref. I17l for the AV18 potential have been also reported for sake of comparison. The results 
obtained by other techniques are also listed. 


Interaction 

Method 

B 

<T> 

vV> 

Pp 

Pd 

Pt =i 

Pt =2 

CD Bonn 2000 

HH (this work) 

26.13 

77.58 

1.454 

0.223 

10.74 

0.0029 

0.0108 


FY [££] 

26.16 

77.59 

- 

0.225 

10.77 

0.0030 

0.0108 

N3LO 

HH (this work) 

25.38 

69.24 

1.516 

0.172 

9.289 

0.0035 

0.0024 


FY [66] 

25.37 

69.20 

- 

0.172 

9.293 

0.0033 

0.0024 


NCSM [20] 

25.36(4) 

- 

- 

- 

- 

- 

- 

AV18 

HH [IT] 

24.210 

97.84 

1.512 

0.347 

13.74 

0.0028 

0.0052 


FY [g] 

24.25 

97.80 

- 

0.35 

13.78 

0.003 

0.005 


FY [69] 

24.223 

97.77 

1.516 

- 

- 

- 

- 


IV. SUMMARY AND CONCLUSIONS 

The A = 3 and 4 nuclear ground states have been studied with non-local two-body potentials using the HH method. 
The variational wave function is written as an expansion over a complete basis, which is constructed in coordinate 
space as a product of hyperradial, HH and spin-isospin functions. The main task is the calculation of the Fourier 
transform of the expansion basis. However, given the properties of the HH functions, this Fourier transform reduces 
to a one-dimensional integral which can be obtained analytically. The application is therefore rather straightforward, 
and the matrix elements of a given two-body interaction can always be reduced to three-dimensional integrals (two- 
dimensional for local potentials). We have presented the results for the CD Bonn 2000 and N3LO NN potential 
models. In both cases we have found very good agreement with the results obtained by other groups. We have also 
given the estimates obtained for various ground state properties and pointed out that the different models predict 
very different (up to a factor 5) percentages of the isospin admixtures in the 3 H and 4 He ground states. 

In recent years, new models of nuclear interaction have been constructed from chiral perturbation theory. They are 
non-local and given in momentum space |38}. The N3LO potential is the first potential of this type fitting the NN 
data set with y 2 /datum ~ 1. Another class of recently developed potentials are the so-called low-q potentials (45l |46j. 
They are “renormalized” two-body potentials V/vat(A;, k'), given in momentum space, where the tail for high values of 
k and k ' has been eliminated. 

It is important, therefore, to have accurate techniques to solve the Schrodinger equation for non-local momentum- 
space potentials. Up to now, only the FY i E3 and the NCSM (20] methods were available for this task. In this 
paper, we have shown that also the HH method can be successfully applied to treat this kind of potentials. 

There are two other important motivations behind this work. The first one is that the HH formalism can be applied 
also to scattering problems. As discussed in Sec. |T| there are still a number of theoretical problems in the N — d, and 
p — 3 He reactions, and it would be very interesting to check whether the use of these new potential models could solve 
these problems. In particular, it appears promising the development of more realistic models of TNIs to be tested 
in A = 3 and 4 nuclear systems. The second motivation is the possibility of the extending the HH method to larger 
systems. The feasibility of such an application would require the solution of several problems, having to do with the 
fast computation of the Raynal-Revai coefficients for A > 4 and the very large degeneracy of the basis. However, 
since these new potentials are softer than the older ones, it should be not too difficult to solve this latter problem. 
Work in both directions is currently underway. 
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